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VARIATIONAL FORMULATION OF RELATIVISTIC 
FLUID THERMODYNAMICS 

L. A. Schmid 

Goddard Space Flight Center, NASA 


ABSTRACT 

As a first step toward generalizing reversible thermodynamics 
from the case of a homogeneous system to that of a system whose 
local velocity may be a function of its position in space-time, a vari- 
ational principle is derived for relativistic reversible adiabatic flow 
of a compressible fluid. This is done by identifying the thermody- 
namic internal energy function for a given sample of the fluid with 
its Hamiltonian function, and then invoking the canonical equations 
of motion. Both in order to bring the rest -mass energy into the for- 
malism, as well as to provide a means of labelling and identifying 
different samples of fluid, it is necessary to introduce a newthermo- 
dynamic variable, which is just the molar initial momentum vector 
of the fluid sample in question. It turns out that this vector is inti- 
mately related to the vorticityof the flow, and if it had been omitted, 
the formalism would have been implicitly limited to a description of 
vorticity-free flow. 

The Lagrangian density, as seen in the fixed laboratory frame, 
that results from identifying the Hamiltonian with the thermodynamic 
internal energy is just the thermodynamic pressure. This must be 
regarded as a function of the generalized coordinates that are canon- 
ical to the particle density, the entropy density, and the initial mo- 
mentum vector (all regarded as generalized momenta). More pre- 
cisely, the pressure is a function of the proper-time derivatives of 
these coordinates. These time derivatives are equal to the molar 
free enthalpy, the rest-temperature, and the initial velocity respec- 
tively. Because, in the laboratory frame, the proper-time derivative 
of a variable is defined as the contraction of the velocity 4-vector 
with the 4-gradient of the variable, the pressure is also a function 
of the fluid velocity. 
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This variational principle yields the correct form of the stress- 
energy tensor for reversible adiabatic flow of a compressible fluid 
(together with the necessary statements of particle and entropy con- 
versation), and automatically gives the expression for the solution 
of Euler's equation of motion for the fluid in terms of the 4-gradients 
of the generalized coordinates. 
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VARIATIONAL FORMULATION OF RELATIVISTIC 
FLUID THERMODYNAMICS 


INTRODUCTION 

The discussion of this paper is in the spirit of well-known attempts 1 to 
bring continuum mechanics within the framework of thermodynamics by treating 
local velocity as just one more thermodynamic variable to be taken into account 
with all the others. The basic approach consists of identifying the appropriate 
thermal energy function with the Hamiltonian of the system, and the corresponding 
canonical equations with the mechanical and thermodynamical equations of motion 
of the system. 

This general approach was first applied to the case of a homogeneous system 
by Helmholtz 2 in 1886, and adapted to Relativity Theory in 1907 by Planck. 3 
Planck's theory was developed before four-dimensional tensor analysis and the 
modern covariance concept had fully evolved. Consequently, although it was form- 
invariant under Lorentz transformations, it fell completely outside the frame- 
work of tensor analysis, which meant that, for all but the simplest applications, 
it was completely unworkable. (Reviews of both the early 4 and recent 5 history 
of Relativistic Thermodynamics are available elsewhere.) 

In 1939 Van Dantzig 6 constructed a manifestly covariant thermodynamics, 
and applied it to fluids, 7 but his work failed to lift the obscurity surrounding the 
intimate three-way relation that binds together thermodynamics, fluid dynamics, 
and the canonical formalism. This relation stems from the fact that, if the right 
choice of variables is made, the thermodynamic energy density function plays the 
role of Hamiltonian density, and the thermodynamic pressure plays the role of 
Lagrangian density. 

The identification of pressure with Lagrangian density had already been 
made in 1908 by Hargreaves 6 for the case of non- relativistic potential flow 
(i.e., vorticity-free flow). Van Dantzig 7 generalized this identification to the 
relativistic case but, although the point was not explicitly made, his proof 
was likewise limited to the case of potential flow because he did not include 
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the variables tnat are necessary lor a completely general description of 
vorticity. (Others have since given relativistic variational principles that 
are free of this limitation, but these principles all involve the imposition 
of constraints, and do not make the identification of the Lagrangian density 
with the thermodynamic pressure.) 

Notation: The analysis will be carried out entirely within the framework of 
Special Relativity. Boldface Latin or Greek letters will designate 4-vectors, and 
light-face characters will designate scalars. An overhead dot will designate 
differentiation with respect to proper-time r, i.e., the time derivative as seen 
by an observer moving with the fluid. Contraction of two 4-vectors will be 
indicated as the dot product of the corresponding boldface characters. Indices 
will be explicitly indicated only in the case of two-index tensors, and when 
indices are indicated, the summation convention will be used. 

Intensive thermodynamic quantities, and extensive quantities that are referred 
to one mole of the fluid, will be designated by capital letters. Thus T and P are 
temperature and pressure respectively, and V, S, U, H, and G are the molar 
volume, entropy, energy, enthalpy, and Gibbs function (free enthalpy) respectively. 
The number of moles per unit volume is n - l/V. Extensive quantities referred 
to unit volume (not unit mass! ) of the fixed laboratory frame will be designated 
by the appropriate lower-case Roman character. For example, u ~ nU is the 
internal energy per unit volume in the laboratory frame. Densities referred to 
the convected fluid frame that is based on coordinate planes embedded in the 
fluid and moving with it will be designated by the corresponding primed letter. 

Thus n' and u' are respectively the molar density and molar energy density 
referred to the convected frame. 


ONE-DIMENSIONAL CANONICAL FORMALISM 

From the point of view of an observer who remains stationary with respect 
to a given sample of fluid and refers all measurements to the converted frame, 
everything can be described as a function of a single var ; " , ble - the proper- 
time t of the sample of fluid under study. Because the fluid appears to remain 
at rest, the fluid velocity v does not enter into such a description. When the 
canonical formalism derived from such an approach is referred to the fixed 
laboratory frame, however, proper-time differentiation must be defined as 
d/d r v • a where 3 is the 4-gradient operator, and this brings v into the formal- 
ism. Thus the development of the one-dimensional canonical formalism 
referred to the convected fluid frame is the first step in arriving at the 
desired variational principle referred to the laboratory frame. 

In reversible adiabatic flow the molar entropy and the total number of 
particles in the fluid are conserved quantities. Our approach will consist of 
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expressing these two conservation laws in terms of two scalar constants of 
motion of the fluid. The internal energy will then be written as a function of 
these two constants of motion and of the proper time. Identifying the internal 
energy with the Hamiltonian of the system and the constants of motion with 
generalized momentum coordinates, we are led to the canonical equations of 
motion. 

In order to arrive at the desired statement of conservation of particles, we 
first note that the molar rest-volume V (not to be confused with the Lorentz- 
contracted molar volume V* = V/l ' where P [ 1 - (v/c) 2 l“ l 2 ) may be written 
V = JV* wherj V' = (V) is the molar volume referred to the convected frame, 
which is a constant of motion, ard is equal to the initial value of V at r = 0, and 
J is the function of r that describes the time-dependence of V that results from 
compression or expansion of the fluid. The Lorentz -contracted molar volume 
is thus V* = V/r = (J/P) V'. Because intervals of laboratory-time dt and proper- 
time d r are related b” dt = T dr we have 


d l) - cd t d V* = c rdr(J/r>dV' r J(cd rdV') 


(1) 


where d U is the element of 4-volume in the laboratory frame and cd r dV* is the 
corresponding 4-volume element in the convected frame. Thus J is just the 
Jacobian of the transformation between laboratory coordinates and converted 
coordinates. 

Using V = JV’, the thermodynamic equation dU = TdS - PdV would become 


dU = TdS - (PJ)dV' - (PV'j)dr 

where S = V* = 0 and U = U(S, V’, r ) would be the thermodynamic potential that 
we could identify with the Hamiltonian. However, because we are dealing with a 
continuum, it is more appropriate to work with densities rather than with molar 
quantities. For this reason, we eliminate V', U, and S in favor of n’, u’, and s’ 
w'here 


1/V' J/V=Jn. 


n' U. 


n'S. 


(3) 
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Making these substitutions in (2), we find 


du' r Tds' + Gdn' - ( P j ) d r 


( 4 ) 


where G = U + PV -TS is the molar Gibbs function (or free enthalpy). Thus 
u' = u' (s', n',r) is a function of two constants of motion and of the proper-time. 

Before identifying u' with the Hamiltonian of the system, we note that (4) 
has two deficiencies which luckily can both be removed by the addition of a single 
term. First, from the relativistic point of view, the rest-mass energy density 
m'c 2 - n'Mc 2 (where M is the molar rest-mass) should not be isolated from all 
other contributions to the energy density. Hence u* should be replaced by the 
total energy density u* n’U that includes the rest-mass energy density. 

The second deficiency of (4) arises from the fact that, if we are to describe 
a fluid, rather than just isolated moles of gas that in no way interact with one 
another, then we must in some way introduce into the formalism parameters 
that label and identify each mole of gas and distinguish it from all others. Be- 
cause these parameters will enter into the formalism, they must have a physical 
significance that is essential to the description of the fluid. Both of these re- 
quirements, labelling and physical significance, are satisfied by the initial 
momentum vector K = (Mv) r = 0 which is the momentum possessed by the mole 
of gas at r = 0. In doing this we are effectively postulating that, the inability 
to distinguish between two or more moles of gas that would result if their K- 
vectors were all equal, represents a physical degeneracy with observable con- 
sequences. (We shall, in fact, see that such a degeneracy corresponds to 
vorticity-free flow.) Thus K , like the molar entropy S, is a preserved fossil of 
the initial conditions of the fluid. The vector K is normalized to the molar mass 
M, i.e., K = (K • K| 1/2 = Me, and so M s K/c can be used as the definition of 
molar mass, and II becomes U U + c(K-K) 12 . There exists an alternative 
procedure for relating U and U that is not only more general, but also closer to 
the spirit of thermodynamical formalism. We may regard U = U(S, n,K ) as the 
basic thermodynamic potential and S, n (or V), and K as the basic variables. We 
then define M as Me 2 (BU/t^K) -K. This definition is consistent with tf = U(S, 
n) + c(K*K ) 1/2 where c(K-K) 1 2 = Me 2 , but it is more general, and is applicable 
regardless of the K -dependence of U. Using this more general definition of M, 
we define the purely thermal energy function U as 


U - U - ( SU/3K ) • K , U - Me 2 . 


( 5 ) 


4 



Although (5) represents the most general way of defining U and M, in this paper 
we shall assume that the K -dependence of 0 is given by c(K*K )* 2 = Me 2 where 
M is a constant parameter. In such a case 5U /aK = ck/K = k /M = v where 
v = (v) t = q is the initial velocity of the mole of gas in question at r = 0. 

If k' * n ’K is the initial momentum density referred to the convected frame, 
the density relation that corresponds to (5) is 


u' * u' - v *k ' = u' - m'c 2 where v e du'/ Bk' = BU/BK . (6) 


Thus we see that the definition of u* in terms of and k ' (or of U in terms of 
& and K) amounts to a Legendre transformation that replaces the variable K ' 
(orK ) with v * B u'/3k' = 3U/BK . 

From (4) and (6) we obtain the basic thermodynamic equation of the fluid: 


du' = Tds' + Gdn' +vdit' -(PJ)dr. 


This is to be compared with the well-known expression for the differential of 
the Hamiltonian E = E(p, q, r): 


dE = 


( ^E/Bp) qT dp + 


(3E/Bq) p T dq + ( BE/Br) p q dr 


= V 1 

L, 


q dp - 



- ( dL /dr) . d r 

q. q 


( 8 ) 


where the Lagrangian L is defined as follows: 


L=L(q,q,r) ■ ^ p( ?E/Bp) q T - E . 

p 


0 ) 
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We now identify u' with E. Note that although iX* = n’t) is a density, it does 
in fact represent the energy of a fixed number of particles, namely the number 
contained in unit volume of the convected frame, and so there is no inconsistency 
in regarding it as the Hamiltonian of a definite dynamical system. 

It turns out that, for consistency, it is necessary to identify the thermo- 
dynamic variables s\ n', and k' with generalized momenta, rather than with 
generalized coordinates. Doing this, and designating the coordinates q that are 
conjugate to the momenta p = (s’, n*, * ' ) by q = (3 , (J ,§ ) respectively, com- 
parison. of (7) and <d) yields the following equations: 


T=3; G = (J; v =5; PJ = (3L/3r) # . < 10 ) 

q. q 


The fact that it' is independent of the coordinates 3, Q, and £ yields the 
desired equations of motion; 


s' = ri' = k ' = 0. which impl ies K = 0. (11) 

the la»t equation resulting from k’ - n'K and n' = k' = 0. 

From (9) we find 


L « 



p(BuVBp) -u' = n'(G + ST + K-v - U) 


= n' (H - U) =n'PV:(n' /n) P = JP 

where use has been made of (3) and (5). Since J dJ/dr = d J/dr , (12) together 
with the last equation of (10), yields (c?P/9r) . = 0, which means that P = P(^' , 

') , 5 ) is a function of the generalized velocities 3 = T, 0 = G, 5 =w alone, and 
not an explicit function of r . For example, in the case of a perfect gas, for which 
P = nRT, where R is the gas constant and y - C p /C v = constant is the ratio of 
the specific heats, the functional form of P is 

P = PoCT/T^^O- 1 )! exp { [G + Mc(v - w) 1/2 ] 'RT } . (13) 


where P„ and T are constants. 

o o 
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Because L = JP, the Lagrangian equations of motion are 


d[3(JP)/Bq] /dr = 3(JP)/3q or d [J( 3P/3q)] /dr = 0 (14) 


where use has been made of the fact that P is independent of the q f s, and J is an 
explicit function of r, being independent of the q’s and q's. To evaluate BP/3q 
we first note that 


P = n(H -U) = n(G + ST + K-w - U) = nG + sT + k *w - u, (15) 

where now the densities are all referred to the laboratory frame. Next we note 
that, from (5) and the relation du = Tds + Gdn, we have 


du = Tds + Gdn + v *dn . : 

Taking the differential of (15) and using (16), we find 

dP = sdT + ndG + k • dv = sd’'' + ndl + k *d§ • (17) 

Using (17) to evaluate BP/?q, we arrive at the following Lagrangian equations 
of motion: 


0 = d( J s ) /dr = s ' ; 0 = d( Jn ) /dr = n' ; 0 = d( Jk )/dr= k' • (18) 


These, of course, agree v/ith the canonical’ equations (11). (If we had identified 
some or ail of the thermodynamic variables with generalized coordinates q, 
rather than with generalized momenta p; this agreement would not have occurred.) 

In the same way that, in arriving at (5) and (6), we noted that the mass density 
m T s n'M could be defined in terms of the k' -dependence of u’, we now note 
that the mass density m = nM (referred now to the laboratory frame) can be 
defined in terms of the -dependence of P: 


me 2 = ( Bu/Bk )*k ~v*k — ( BP -'Bg ) • 5 • 


(19) 
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Thus, the definition of the molar mass M may be taken to be 


M = n" 1 (BP/?§ ) -§/c 2 = (dP/Bf))" 1 /c 2 . < 20 ) 

VARIATIONAL PRINCIPLE FOR FLUID 

The Lagrangian equations (14) can be obtained from the following variational 
principle: 


0 = 



JPdr. 


This refers to the fluid contained in unit volume of the convected frame. If the 
integrand were JPd r dV\ the principle would refer to the sample of fluid con- 
tained in the volume dV\ Since the fluid contained in each volume element dV' 
must individually and independently satisfy the requirement 


0 = 5 


T 



JPdrdV' 


then it must follow that 


0 = 8 f JPdrdV' 

J r.v' 

where now the integration extends over V' as well as over r. Thus, referring to 
(1), we arrive at the following variational principle for the fluid: 


0 = 5 


i PJcdrdV' 

Lv 



( 21 ) 
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where d 0 = cdtdxdydz is the 4-volume element in the laboratory frame, and 0 
is the 4-volume occupied by the fluid between the specified initial and final times. 

Because the proper-time r is no longer the independent variable, the opera- 
tion d/dr must be defined in terms of the fluid velocity v as d/d r = vd where 
d is the 4-gradient operator. Because v must remain normalized (v*v = c 2 ) 
during the variation process, it must be parameterized in some way so that the 
normalization will be guaranteed. This can be done most conveniently by intro- 
ducing a vector p whose direction is v : 


v = cp /p where p T (p‘p) 1/2 - (22) 

It will turn out that the norm p does not appear in any of the Euler-Lagrange 
equations resulting from (21). The components of p are to be regarded as 
generalized coordinates, rather than as velocities. Using the definition (22) 
for v we have 


G = (} = v-3Q = p -1 c p • • 


Similar expressions define 3 and 5 . 


(23) 


The more detailed statement of the variational principle given in (21) is 


0 = 



(5q)[(BP/Bq)- 3-p]dU + 



( Sq)p • 


(24) 


where p - BP/B (aq) and is an element of the hyper-surface $ that bounds 
the 4-volume l) over which the integration is carried out. The variational principle 
is thus equivalent to the requirement that the Euler-Lagrange equations 3 *p = 

BP/B q be satisfied, and that the variables have definitely assigned values on the 
boundary so that Sq = 0 on k . 

Referring to (17) and (23), the calculation of the generalized momenta can be 
illustrated by the case forQ : 

P n - BP/B(dQ) = (BP/BG) [BG/B(3Q>] = nv . (25) 
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Similarly we find p n= Snv, p^ k = nK k v* , and p*, k =0. Thus the Euler- Lagrange 
equations corresponding to variation of (J , 3 , and g respectively are 


0 r d • ( vn) = 3 ■ ( vnS ) = 3 • ( vnK ) . 


.'26) 


Variation of p yields the following equation: 


0 = 3P/3p= ^ ( 3p /^)[3(v3q)/ 3 p] 

q 


= (c 



( BP/Bq) [3q - v (q/c 2 )] 


= (cn /p) [3Q + S33 + (3g) -K - (M + H/c 2 )v] (27) 


or 


(M + H/c 2 ) v = 3Q + S33 + (3§)-K. < 28 ) 

Equations (26) are just the required conservation laws. Using the first, the 
second and third could also be written as S = K =0. Equation (28) is effectively 
the formal solution of the fluid equation of motion (Euler's equation), and thus 
amounts to a statement of conservation of energy-momentum. That this is true 
can be verified by evaluating the stress-energy tensor w which, since our 
Lagrangian density is P, is given by 


w 


i 

k 


^ P j 3 k q - PSj[ = ^ [3P/3(3.q)] 3 k q - psj 

q q 


(29) 


which, if the Euler- Lagrange equations 3 • p = 3P/3q are satisfied, automatically 
satisfies the equation 


3jwi =-(3P/3x k ) qii =0. 


(30) 



where we have postulated that the pressure function possesses no explicit de- 
pendence on the space-time coordinates. Using the expressions for the p's that 
were given following (25), and making use of (28), we find that (29) becomes 

w l =nvi [3 k Q + S5 k 3 + (\S)*K]-P8i 

= n(M + H/c 2 )v>v k - Pgj • (31) 

Making use of the first equation of (26), equation (30) becomes 

d [(M + H^c 2 ) v]/dr= n~ 1 #P. (32) 

This is just Euler's equation for the fluid, and may be regarded as the determining 
equation for the molar energy-momentum vector (M + H/c 2 ) v . But (28) gives an 
explicit expression for this vector in terms of the 4-gradients of the canonical 
coordinates, so, as previously remarked, (28) constitutes the formal solution of 
Euler's equation. 

It should be noted, incidentally, that when K becomes constant over any 
region, the term (B§ )*K in (28) becomes the gradient of a scalar, and this cor- 
responds to vorticity-free flow 9 in this region. As previously noted, this 
physically observable effect is characterized by a degeneracy resulting from the 
fact that the labelling vector K is indistinguishable for neighboring samples of 
fluid. If the vector K, and hence § , had never been introduced into the formalism, 
and we had instead introduced the rest-mass energy Me 2 simply by replacing 
0 by % where d^/dr = G + Me 2 , we would have arrived at a variational principle 
implicitly restricted to the case of vorticity-free flow. 
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